Tunneling anisotropic magnetoresistance driven by resonant surface states: 
First-principles calculations of a Fe(OOl) surface 
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Fully-relativistic first-principles calculations of the Fe(00f) surface demonstrate that resonant 
surface (interface) states may produce sizeable tunneling anisotropic magnetoresistance in magnetic 
tunnel junctions with a single magnetic electrode. The effect is driven by the spin-orbit coupling. It 
shifts the resonant surface band via the Rashba effect when the magnetization direction changes. We 
find that spin-flip scattering at the interface is controlled not only by the strength of the spin-orbit 
coupling, but depends strongly on the intrinsic width of the resonant surface states. 

PACS numbers: 72.25.Mk, 73.23.-b, 73.40.Gk, 73.40.Rw 



Spin-dependent properties of magnetic surfaces and in- 
terfaces have long been a subject of vigorous research. 
Recent interest is triggered by the advent of magnetoelec- 
tronics, a technology aimed at harnessing the electrons 
spin in data storage and processing. Typical devices 
utilize heterostructures composed of magnetic and non- 
magnetic materials [l[ . One approach recently suggested 
takes advantage of spin-orbit coupling (SOC), where the 
resistance of a nanostructure is controlled by rotating its 
magnetization with respect to the interface. For a tun- 
nel junction with only one ferromagnetic electrode this 
effect is called tunneling anisotropic magnetoresistance 
(TAMR) [1 0,0,0,11. It is different from the usual tun- 
neling magnetoresistance (TMR) observed in magnetic 
tunnel junctions 0] with two ferromagnetic electrodes 
when their orientation is switched between the paral- 
lel and antiparallel states. It is also different from the 
usual anisotropic magnetoresistance (AMR) typical for 
bulk ferromagnets [|| and from ballistic anisotropic mag- 
netoresistance (BAMR) which may occur in ferromag- 
netic nanocontacts [9(. The TAMR occurs in the tunnel- 
ing transport regime because SOC makes the electronic 
structure anisotropic. Practical advantages of TAMR 
over AMR are: (1) the tunneling process filters out a 
fraction of the electronic phase space, and (2) SOC is 
usually stronger at interfaces. Both these features result 
in a larger effect. An advantage of a TAMR device over 
a conventional magnetic tunnel junction is that only one 
magnetic interface is necessary. 

Experimental demonstrations of TAMR to date have 
employed diluted magnetic semiconductors (GaAs:Mn) 
as magnetic electrodes @, 0, 0, 0, 0). The TAMR effect 
is due to the anisotropy of the tunneling density of states 
in the valence band of GaAs induced by SOC 0> 0]- It 
was suggested that metallic alloys, such as CoPt, with 
large magnetocrystalline anisotropy may be used as elec- 
trodes in TAMR devices [13]. TAMR was also observed 
by scanning tunneling spectroscopy (STS) of thin Fe films 
on W(110), as the result of a group velocity change from 



the band splitting induced by SOC [11], 

In this Letter, we propose a different route to achieve 
large TAMR with a metallic electrode. It is well known 
that many transition-metal surfaces, as well as their in- 
terfaces with insulators, exhibit electronic bands that are 
localized at the surface or interface. If such a band 
does not mix with bulk states it is called a surface 
band. If it mixes weakly with bulk bands, it broad- 
ens and becomes a resonant surface band. Interface 
bands usuall y c ontribute strongly to the tunneling cur- 
rent [H, El tLI EE El- On the other hand, when the 
magnetization is rotated, these bands shift in the inter- 
face Brillouin zone due to the Rashba effect [13] pro- 
duced by SOC at the interface. For non-magnetic sur- 
faces the Rashba effect lifts the spin degeneracy of the 
surface states Q; for magnetic surfaces it leads to an 
asymmetric shift of the non-degenerate band [l^ . The 
idea pursued in this Letter is to use this Rashba shift in 
the TAMR device. To illustrate this approach, we choose 
the Fe(001) surface, because it is one of the most studied 
transition-metal surfaces supporting a surface minority- 
spin resonant band at the Fermi level [H ■ We show that 
the change of the magnetization direction both in and 
out of the surface plane results in a sizable change in 
the tunneling conductance, which may be observed using 
STS. We also identify conditions controlling the strength 
of spin-flip scattering at the interface. 

To study TAMR produced by electron tunneling from 
the Fe(001) surface through vacuum we need a counter- 
electrode for closing the electric circuit. For this pur- 
pose it is convenient to use a metal which does not fil- 
ter the electrons by the transverse wavevector k||. A 
non-magnetic bec Cu electrode has a spin-independent 
free-electron-like band structure and a featureless surface 
transmission function [l3j |. which makes it an ideal spin 
detector, similar to the STS tip in the Tersoff-Hamann 
theory [2fjj. The semi-infinite Fe and Cu leads are sep- 
arated by approximately 1 nm of vacuum (6 monolayers 
of empty spheres) . The conductance is calculated using 



the Landauer-Biittiker approach plj implemented within 
the fully-relativistic tight-binding linear muffin-tin or- 
bital (TB-LMTO) method [13]. Charge self-consistency 
is achieved using scalar-relativistic TB-LMTO calcula- 
tions for Fe and Cu surfaces treated using supercells with 
12 metallic monolayers. 

The numerical technique is based on the Green's func- 
tion representation of the (TB-LMTO) method in the 
atomic spheres approximation j23|. Within the relativis- 
tic formulation of the local spin density approximation 
in which only the spin component of the current density 
is taken into account [241 ] . inside each atomic sphere we 
solve the Kohn-Sham Dirac equation [25[. The technique 



is similar to Refs. 1261 . 1271 : the primary difference here 
is that we use third-order potential functions 28]. The 
Green's function of the layered system is constructed by 
a principal-layer technique |29l|. and the conductance is 
calculated similar to Ref. [3fJ. The surface Green's func- 
tions of the leads are constructed scalar-relativistically, 
and therefore the conductance G is represented as a sum 
of four spin components G aa i = (e 2 /h)T aa i (3lj |. where 
Tfjfj/ is the transmission function integrated over kii (a 
uniform 200x200 mesh was used for this integration). 

Figs, nja-c) show T aa i for three magnetization di- 
rections. The energy dependence represents the linear- 
response conductance in the rigid-band model, which ap- 
proximately reflects the effects of alloying. It is seen that 
Tfj exhibits featureless free-electron-like energy depen- 
dence. However, is nonmonotonic and dominates in 
the energy range between —125 meV and 25 meV. The 
TAMR ratios (T fi - li 100 !)/^ 100 ! are shown in Fig. [Hi 
for both out-of-plane (n = [001]) and in-plane (n = [110]) 
magnetization orientations as a function of bias voltage 
[32I ]. In both cases the TAMR has a spectacular change 
of sign close to the Fermi level and reaches ±15-20% at 
the bias voltage of =f50 mV. 

To explain the origin of large TAMR effect, we fo- 
cus on energy —15 meV where the conductance has the 
strongest spin asymmetry. Fig. [2] shows spin- and k||- 
resolved surface densities of states (DOS) and tunneling 
transmission functions. The left six panels correspond 
to the out-of-plane magnetization (n — [001]), and the 
right half to the in-plane magnetization (n = [100]). The 
resonant surface band is responsible for the bright four- 
petal-flower features in the minority-spin surface DOS 
( -/Vj 1 and Nf ) and in the minority-spin transmission (Tj| 

and t\). This band is dominated by the minority-spin 
surface states, which mix weakly with bulk bands. A 
central point is that SOC can strongly enhance this mix- 
ing, in both spin-diagonal and spin-mixing components. 
In particular, consider the surface state lying on the T X 
line with k|| = (k x ,0). In the absence of SOC , these 
eigenstates have definite parity with respect to reflection 
in the y = plane. The surface band is even, while 
the minority-spin bulk band is odd. By symmetry these 
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FIG. 1: (a-c) Spin-resolved integrated transmission T CT(T ' for 
the Fe(001) surface with a Cu counterelectrode as a function 
of energy. Magnetization is along (a) [100], (b) [110], and (c) 
[001] directions. The Fermi level is at zero energy, (d) In- 
plane and out-of-plane TAMR as a function of bias voltage. 



states cannot mix, and the surface state remains local- 
ized. The SOC does not conserve this parity and mixes 
the surface state with both minority-spin and majority- 
spin bulk states. The surface state is thus transformed 
into a surface resonance. In our case this occurs at spe- 
cial k|| points. In general, if for a given spin the surface 
band lies within a gap of bulk bands, the SOC converts 
this localized surface band into a resonant band. 

As is evident from Fig. [3J the surface bands depend 
on spin orientation (compare to ivj' and to N')). 
This can be attributed to the Rashba effect [171 ] . The 
effective spin-orbit shift of electron energy is given by 



Ae(k||) — a (z x ky) • s, 



(1) 



where z is the unit vector normal to the surface. The 
electron spin s is aligned with the magnetization M. For 
M || z, the Rashba shift is zero throughout the surface 
Brillouin zone, while for M || x it is positive for k y > 
and negative for k y < 0. The asymmetric shift of the 
surface bands is reflected in the loss of fourfold symmetry 
in the right six panels of Fig. [2] and in particular in 
the loss of mirror symmetry in y. A similar effect was 
discussed for the Gd(0001) surface [3]. 

Comparison of to N^-, as well as 7V|' to ivj , 
in Fig. [2] indicates that the admixture of majority-spin 
states to the surface band is of the order of 1% (note the 
difference in scales). However, the spin- flip components 
of the transmission are quite pronounced. For example, 
in certain areas of the Brillouin zone both T-\ ^ and 
are comparable to the majority-spin component . For 
M || z, a portion of the resonant surface band lies close to 
the r point and adds a large contribution to the minority- 
spin conductance. When the magnetization is rotated to 
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FIG. 2: Spin components of the DOS at the Fe(001) surface monolayer and the transmission functions T aa i for the Fe(001) 
surface with a Cu counterelectrode at 15 meV below Ef- Figures are resolved by kn with abscissa along [100] and ordinate 
along [010]. The left six panels correspond to the magnetization normal to the surface (labeled by superscript _L), and the right 
half to the in-plane magnetization aligned along [100] (labeled by superscript ||). Some panels are given in a logarithmic scale. 



M || x these states shift, and the conductance is reduced. 
This is the origin of the large TAMR effect seen in Fig. [I] 

The spin-flip components of the transmission function 
Tjj and Tj_j shown in Fig. [T] display a nonmonotonic en- 
ergy dependence and are generally quite small compared 
to the spin-conserving components. Surprisingly for the 
in-plane magnetization, has a pronounced maximum 
at Ep — 0.1 eV, just above the bottom of the resonant 
surface band, which extends to higher energies (Fig. [T^,). 
Here the spin-flip and spin-conserving contributions are 
comparable. Notably the peak appears only for the in- 
plane magnetization, and the spin-flip process is strongly 
asymmetric: 3> T^. Fig. [3] shows ky-resolved spin- 
flip transmission function Tj j,(kii ) for both magnetization 
orientations. The resonant surface bands are seen as four 
small ellipses along the T X lines. All the difference in 
for the two orientations accrues from these four ellipses, 
which clearly indicates that the large spin-flip conduc- 
tance is entirely due to the resonant surface states. 

To elucidate the origin of the resonant spin-flip trans- 
mission we consider a simple tight-binding model. We 
assume that there is a single majority band, with sur- 
face Green's function GS* in the absence of spin-orbit 
coupling, and a minority surface band -E s (kii), broad- 
ened by a hybridization with the bulk minority band. 
The broadening is included through parameter 70 that 
changes the corresponding Green's function into G? 1 = 

(_E — E s (k\\) + 170) . The spin-orbit interaction Vso 
mixes spin channels resulting in the "dressed" surface 
Green's function G a a' ■ 



We apply the Landauer-Biittiker formalism [2l| to 
the case of two coupled spin channels, treating the sur- 
face monolayer as a "conductor" and everything else as 
"leads." We assume that the tunneling probability is 
small and that the second electrode is non-magnetic. 



Then 



E 

<j<j' 



T ™' =4^ImS CT |G ffff ,| 2 ImE^ (2) 



where E^ is the spin-dependent self-energy of the surface 
layer due to interaction with the magnetic electrode, and 
E„ the spin-independent self-energy due to interaction 
with the nonmagnetic electrode through vacuum. In- 
cluding Vso through the Dyson equation, we find the 
spin-flip component of the Green's function 



G tl = 



V so G° t 



11 



£-£ s (k||)+z7o-|l/so| 2 GO T - 



(3) 



Using the identity ImGi^lG^I 2 ImS-p, and writing 
— |Vso| 2 ImG|| as 7, we find for the spin-flip transmis- 
sion coefficient: 



T Tl( k ||) K 



[E - E s (k|| ) - A] " + (70 + 7)2 



(4) 



where A = | Vso 1 2 ReGijL is the shift of the surface band. 
Finally, we add an energy delta-function and integrate 



over k|| [33j to obtain 



T u {E)kN.{E) 



7 



7o+7' 



(5) 



where N S (E) is the DOS of the surface band. 

The parameters 7 and 70 can be interpreted as broad- 
ening of the minority surface band due to the coupling 
to the majority band (through Vso) and to the minority 
band (through both Vso and hopping matrix elements), 
respectively. Eq. (O shows that the spin-flip transmission 
depends on 7/70, and is large when 7/70 is large. As 
was shown above, along the T X lines there is no mixing 
at all with the minority band when Vso = 0. Near the 



4 





FIG. 3: Spin-flip transmission Tf± for magnetization along (a) 
[100] (in-plane), and (b) [001] (out-of-plane), at -0.102 meV. 



bottom of the interface band, the surface band shrinks 
to four pockets along these lines (Fig. [3|) , and hence 70 
is small for all the surface states at this energy. 

The directional dependence of (Fig. [3]) comes from 
the relative magnitude of the spin-orbit contributions to 
7 and 7q. For the perpendicular orientation, M || z, V$o 
mixes the interface states primarily to minority-spin bulk 
states, contributing mainly to 70. Therefore 7 <S 70, and 
the spin-flip transmission is small. For M || x, spin-orbit 
contributions to 70 and 7 are of the same order, and 
resonant spin-flip transmission sets in. As follows from 
the Dyson equation, the second spin-flip component 
is smaller than Tfj by a factor of order 70/Wf, where Wj 
is the bandwidth of the majority bulk band. Thus, the 
model predicts <C in perfect agreement with the 
detailed calculations, as seen in Fig. [T] 

To conclude, we have investigated the electronic struc- 
ture and tunneling from the Fe(001) surface to demon- 
strate that large TAMR may be achieved by utilizing the 
Rashba shift of interface resonant states. TAMR values 
of up to ±20% were predicted for the Fe(001) surface 
at small bias voltages; this conclusion may be checked 
experimentally using STS with a non-magnetic tip. We 
also found a spectacular resonant spin-flip transmission 
near the interface band edge, and showed that spin-flip 
scattering at the interface depends strongly on the in- 
trinsic broadening of the resonant band. The results are 
generally applicable to magnetic surfaces and interfaces 
carrying interface states and suggest that large TAMR 
may be achieved using metallic electrodes. 

The work at UNL was supported by the Nebraska Re- 
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